NEW  YORK  UNIVERSfTY 
AFCRC-TN-57-209  INSTITUTE  OF  MATHEMATiC  • 

ASTIA  DOCUMENT  No.  AD  133673  ,^  ...       ,    i:'^^^^^ 

;25  Wavarly  Place,  tiew  york  3,  N.  Y. 


■aEET  PR^x 
^>  -    J^^      NEW   YORK   UNIVERSITY 

^   I   In   IJL    ^     Institute  of  Mathematical  Sciences 


\i^     ^^      ^      Division  of  Electromognetic  Research 
•^  C  CCX  A 


RESEARCH    REPORT    No.    CX-31 


Excitation  of  Molecular  Vibration  by  Collision 


MORRIS  SALKOFF  and  ERNEST  BAUER 


CONTRACT  No.  AF  19(604)1705 
JULY,    1957 


NEW  YORK  UNIVERSm 

Institute  of  ^&theInatical  Sciences 
Division  of  Electromagnetic  Research 

Research  Report  No.  CX-31 


EXCITATICN  (F  MOLECUIAR  VIBRATICK  BY  COLLISICN 
Morris  Salkof f  and  Ernest  Bauer 


Morris  Salkoff 
Ernest  Bauer 


Morpis  Kline 
Project  Director 


July,  19^7 


The  research  reported  in  this  document  has  been  sponsored 
by  the  Geophysics  Research  Directorate  of  the  Air  Force 
Cairibridge  Research  Center,  Air  Research  and  Development 
Command,  under  Contract  No.  AF  19(60U)1705« 


-  1  - 
Abstract 

We  have  calculated  the  cross  section  for  excitation  of  the  first  vibra- 
tional level  in  the  collision  of  two  hydrogen  molecules  in  the  ground  state.  A 
perturbation  scheme  must  be  used  because  of  the  many-body  problem  involved.  The 
cross  section  is  calculated  by  a  partial  wave  analysis.  V/e  consider  the  collision 
of  two  Hp  molecules  because  the  intermolecular  potential  is  reasonably  well  known 
for  this  case. 

Because  of  the  large  masses  of  the  particles,  the  de  Broglie  wavelength 
of  the  system  is  much  smaller  than  the  range  of  the  intermolecular  potential.  Con- 
sequently many  values  of  the  angular  momentum  must  be  considered.  We  find  that  the 
major  contribution  to  the  total  cross  section  comes  from  large  values  of  the  angular 
momentum  -  ^  of  the  order  of  10  to  20  -  and  that  values  of  angular  momentum  ixp  to 
about  80  must  be  considered.  In  order  to  solve  the  Schrodinger  equation  for  arbi- 
trary values  of  angular  momentum  with  a  Morse  potential  form  for  the  intermolecular 
potential,  the  high  speed  computer  Univac  was  used. 

In  the  energy  range  of  interest,  0.5  to  2  ev,  the  total  cross  section 
is  obtained.  At  energies  greater  than  5  ev  the  values  for  the  cross  section  are 
larger  than  the  maximum  permissible  if  particle  flux  is  to  be  conserved,  that  is, 
our  pertxirbation  scheme  breaks  down  at  these  high  energies.  However,  the  cross 
section  for  energies  greater  tJian  5  ev  is  not  too  significant  since  electronic 
transitions,  which  we  have  neglected,  become  appreciable. 

We  find  the  value  of  the  cross  section  to  be  very  sensitive  to  the  pre- 
cise values  of  the  potentials  used.  This  result  is  to  be  expected  since  the  classi- 
cal turning  points  of  the  motion  occur  at  val.ues  of  internuclear  separation  for 
which  the  potentials  are  steeply  varying  functions  of  nuclear  separation.  Conse- 
quently the  calculation  of  an  inelastic  scattering  cross  section  for  most  molecular 
sj^tems  -  for  wiilch  the  details  of  the  intermolecular  potential  are  not  known  - 
must  be  considered  approximate,  at  best. 
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!•  Review  of  previoias  work 

1.1  Description  of  problem 

1.11  Introduction 

The  characteristic  collision  processes  in  a  molecular  gas  in  the  tempera- 
ture range  from  100  °K  to  5000  °K  are  elastic  scattering,  and  inelastic  scatter- 
ing which  involves  changes  in  vibrational  and  rotational  quantvun  numbers  only. 
Electronic  transitions  diiring  collision  will  take  place  to  a  significant  extent 
only  at  much  higher  temperatures,  E  >5  ev( corresponding  to  60,000  K)« 

For  most  molecules,  the  spacing  of  the  rotational  energy  levels  is  small 
(of  the  order  of  .05  ev),  so  that  changes  in  internal  rotational  states  on  colli- 
sion are  very  likely  at  room  temperature  (300  K  corresponds  to  .03  ev)» 

We  are  concerned,  however,  with  the  excitation  cf  vibrational  energy 
states  on  collision.  The  spacing  of  the  vibrational  states  is  of  the  ordsr  0.1  to 
0.3  ev  for  most  molecules.   Since  the  inelastic  scattering  cross  section  near  the 
threshold  depends  exponentially  •- J  on  the  energy  transferred  from  translational  to 
vibrational  degrees  of  freedom,  the  cross  section  will  be  very  small  near  the  thres- 
hold energy  and  increase  rapidly  tdth  energy  above  threshold. 

At  temperatures  below  5000  K, electronic  transitions  do  not  occur  to  any 
appreciable  extent,  and  so  the  Bom-Oppenheimer  *- -*  approximation  is  very  good.  In 
this  approximation  the  nuclei  of  the  molecules  are  assumed  to  move  so  much  more 
slowly  than  the  orbital  electrons  that  the  collision  is  'adiabatic'  in  the  sense 
that  the  electrons  have  time  to  adjust  themselves  without  transitions  to  the  change 
in  potential  energy  during  the  collision.  The  nuclei  'see'  only  the  average  field 

produced  by  the  motion  of  the  electrons.  The  Born-Oppenheimer  approximation  is  good 

■a. 

e.g.,  for  Np,  the  spacing  is  0.29  ev,  for  Op,  0,2  ev.  For  hydrogen,  which  is  of 

coTirse  somewhat  special,  the  energy  difference  between  the  ground  and  the  first 
excited  states  is  0.51;  ev. 
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so  long  as  the  nuclei  do  not  have  sufficient  energy  to  induce  electronic  transitions. 
The  intermolecular  potential  can  then  be  defined  as  a  function  of  nuclear  separation 
only.   The  potential  energy  at  each  value  of  intemuclear  separation  is  calculated 
assuming  stationary  electronic  motion. 

Inelastic  collision  cross  sections  are  needed  to  calc^jlate  rates  of 
chemical  reactions  and  of  vibrational  relaxation.  The  latter  mechanism  explains 
the  dispersion  of  high  frequency  so\md  in  di-  and  polyatomic  gases.  The  inelastic 
scattering  cross  section  is  also  needed  in  the  calculation  of  three-body  recom- 
bination coefficients'--',  and  the  accomodation  coefficient  of  a  gas'--'. 

1.12  Formulation 

Call  the  first  moleciile  AB,  the  second  C,  Tjie  second  molecule  is  con- 
sidered to  be  a  simple  particle  of  appropriate  mass.  The  motion  of  the  center  of 
mass  is  eliminated  from  the  Schrodinger  equation  for  the  system  by  using  the  fol- 
lowing coordinates: 

(1)  P  =  r^ 

t  =  (MaV  ^V  "c^c^/^a"  V  ^c  • 

These  coordinates  are  shown  in  Figure  1.  Here,  p  is  the  internuclear  separation 
of  the  molecule  AB  and  describes  the  internal  vibrational  motion;  r  is  the  distance 

j5 r  ^ 

A  •   S  6 

^ p ^ 


Figure  1 

from  the  second  molecule  C  to  the  center  of  mass  of  AB.  This  describes  the  relative 
translational  motion  of  AB  and  C.  R  is  the  coordinate  of  the  center  of  mass  of  the 
system  as  a  whole  and  need  not  be  considered  explicitly. 
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The  reduced  masses  then  are 

The  coordinates  (l)  are  particularly  suited  to  the  present  probleirit  It 
is  easily  shown  that  the  kinetic  energy  of  the  whole  system  reduces  to  a  s-ura  of 
squares  in  this  coordinate  system  ^  A  This  is  true  whether  the  kinetic  energy  be 
calculated  classically  or  quantum  mechanically. 

There  are  now  two  major  difficulties  in  the  formulation  of  the  collision 
problem  for  inelastic  scatterings  (a)  since  we  are  dealing  with  a  many-body  pro- 
blem, some  suitable  perturbation  scheme  must  be  devised^  and  (b)  an  approximation 
to  the  intermolecular  potential  must  be  found.  The  perturbation  scheme  we  use  con- 
sists in  separating  the  intermolecular  potential  into  two  parts:  an  unperturbed 

O  .A        Jk  t 

part  V  which  is  additive  in  the  two  coordinates  p  and  r,  and  a  perturbation  V 
which  couples  the  internal  (vibrational)  and  external  (translational)  coordinates 
and  thus  induces  the  transition,  '^he   perturbation  potential  is  assumed  to  be  much 
smaller  then  V  .  Thus  we  write 

(3)  V  =  V°+  v' 

ih)  V°  =  V?(?)  +  V°(?) 

J.         o 

(5)  v'  =  v'(p,  ?)  . 

Here  V.  is  the  internal  potential  energy  of  the  molecule  AB,  and  Y°  is 
the  external  (translational)  potential  energy  between  the  two  molecules,  V  is  the 

perturbation  potential.  It  is  clear  that  this  perturbation  approach  will  be  valid 

'  o 

only  so  long  as  V  is  actually  much  smaller  than  V  .  This  will  be  true  only  if  the 

molecules  do  not  approach  each  other  too  closely  during  the  collision.  This  con- 
dition is  violated  at  high  energies,  and  then  the  perturbation  scheme  breaks  down. 
This  result  is  discussed  in  Section  2, 
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The  interraolecvilar  potential  V  is  not  well  known  quantitatively  except 
for  the  simplest  system,  such  as  H--  Hp.  For  this  reason  some  approximation  to 
the  intermolecular  potential  must  be  made.  The  general  features  of  the  inter- 
molecular  potential  are  well  known.  At  large  values  of  r  the  potential  is  es- 
sentially the  attractive  van  der  Waal's  potential  which  varies  as  r  .  At  small 
distances  the  Coulombic  repulsion  of  the  electrons  is  dominant  and  there  is  a 
steep  repulsive  barrier.  At  intermediate  distances  an  exponential  form'--'  will 
be  more  correct  than  a  power  series.  Although  the  general  features  are  well  known, 
it  is  not  possible  to  calculate  analytically  the  intermoleciilar  potential  for  any 
but  the  simplest  systems.   The  angular  dependence  of  the  intermolecular  potential 
has  been  calculated  for  IL-  Hp,  but  not  for  any  more  complicated  molecular  systems. 
Because  of  these  difficulties,  one  attempts  to  fit  a  spherically  symmetrical  poten- 
tial with  two  or  more  adjustable  parameters  to  available  data  on  the  molecular 
system.  This  work  will  be  discussed  in  Section  1.3* 

The  major  mathematical  difficulty  arises  in  the  solution  of  the  Schrodinger 
equation  for  even  the  simplest  form  for  the  intermolecular  potential.   The  radial 
equation  will  involve  the  centrifugal  potential  term +i  i  ii+   l)/2M  r   and  is 
consequently  intractable  for  even  simple  forms  for  the  potential  V  .  Because  of 
the  large  masses  of  the  particles  and  the  relatively  long-range  potential,  many 
values  of  angular  momentum  J?  must  be  considered,  and  this  considerably  complicates 
the  calculation.  It  is  usual  to  perform  the  calculation  only  for  the  case  /  =  0, 
but  past  work  has  shown  the  inadequacy  of  this  approximation'-  -^ ,     This  is  further 
borne  out  in  the  present  work. 
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1,2  Previous  calculations  of  inelastic  cross  sections 
1.21  Work  of  Zener,  Jackson  and  Mott. 

An  important  paper  on  the  excitation  of  molecular  vibration  by  collision 
was  written  by  C«  Zener*-  •^.  The  techniques  developed  there  for  approximating  the 
complicated  many-body  problem  have  been  the  starting  point  of  most  later  calculations. 

Zener  approximates  the  intermolecular  potential  by  the  exponential  form 
e   .  This  is  the  asymptotic  form  of  the  resonance  energy  between  tw-o  symmetrical 
diatoraic  molecules  where  x  is  the  separation  between  the  centers  of  gravity  of  the 
two  molecules.  This  is  of  course  a  highly  simplified  form  for  the  intermolecular 
potential  wMch  ignores  the  attractive  part  of  the  potential.  At  the  energies  of 
interest,  however,  it  is  a  sufficiently  good  approximation  and  moreover  is  rela- 
tively tractable  mathematically. 

Only  collinear  collisions  can  be  considered  since  not  even  tJie  simple 
exponential  approximation  to  the  intermolecular  potential  yields  a  Schrodinger 
radial  equation  that  is  soluble  for  arbitrary  angular  momentum  i  •  This  approxi- 
mation has  been  used  by  many  writers  since  Zener  although  it  is  easily  shown  that 
for  most  systems  of  interest,  many  values  of  Jl    will  contribute  to  the  inelastic 
cross  section, 

A  brief  consideration  of  the  problem  from  the  viewpoint  of  classical 
mechanics  gives  the  folloviing  qualitative  picture  of  the  collision j  when  the  time 
of  collision  is  greater  than  the  period  of  oscillation  of  the  molecule,  the  col- 
lision will  be  'adiabatic'  That  is,  the  colliding  atom  or  moleciile  will  tend  to 
act  only  on  the  center  of  gravity  of  the  struck  molecule,  and  hence  the  probability 
of  a  transition  to  an  excited  ^/Ibrational  state  will  be  small.  Zener  calcxilates 
the  ratio  of  the  time  of  collision  to  the  period  of  vibration  using  a  simplified 
potential  and  classical  mechanics.  The  result  shows  that  for  most  systems  this 
ratio  is  large,  and  thus  we  can  expect  a  small  inelastic  cross  section* 
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These  considerations  suggest  the  use  of  a  perturbation  scheme  to  calcu- 
late the  probability  of  transition  to  an  excited  vibrational  state.  An  approxima- 
tion to  the  exponential  potential  is  used  in  order  to  solve  the  radial  equation 
for  i  =  0  in  terms  of  simple  functions.  Zener  calculates  the  probability  of 
deactivation  of  Np  from  the  first  excited  vibrational  state  to  the  ground  state 
on  collision  with  He  at  room  temperature.  The  result  is  6  •  10  • 

Jackson  and  Mott'--^J  perform  a  calciilation  very  similar  to  that  of  Zener 's 
and  apply  the  results  to  the  calculation  of  the  accomodation  coefficient  of  a  gas« 
The  calculation  proceeds  in  exactly  the  same  way  and  gives  the  same  result  for  the 
transition  probability.  However,  the  authors  are  able  to  solve  the  radial  equation 
for  1=0  with  a  simple  exponential  form  for  the  intermolecular  potential,  in 
terms  of  Bessel  functions  of  imaginary  order.  This  makes  it  possible  to  write  down 
the  transition  probability  in  terms  of  simple  functions  for  this  choice  of  potential. 

Devonshire  "-^-l  and  Zener  ^-^  carried  out  a  slightly  different  calculation 
which  takes  into  account  the  attractive  part  of  the  intermolecular  potential.  ^Iiey 
were  able  to  calculate  the  transition  probability  for  an  intermolecular  Morse  po- 
tential of  the  form 
(6)  V  .  .D[2e-(-^V2a.3-(r-b)/a]  ^ 

TMs  potential  is  zero  at  r  =  b  -'2a(ln  2)  and  has  a  mimimura  at  r  =  b. 
It  is  physically  more  realistic  than  a  single  exponential  even  though  it  does  not 
behave  correctly  for  large  r.  At  large  r  we  expect  the  potential  to  go  to  zero  as 
r  .  '-Che  result  obtained  with  this  potential  does  reduce  to  that  of  Jackson  and 
Mott  in  the  limit  D  approaching  zero. 

Because  the  radial  equation  can  not  be  solved  in  closed  form  for  the 
potential  (6)  even  for  /(  =  0,  a  semi -classical  method  must  be  used  in  order  to 
obtain  the  transition  matrix  element  of  the  perturbation  potential.  This  method 
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agrees  with  a  q-uantiun  mechanical  restilt  only  in  the  limit  of  high  energies.  V/e 
shall  see  in  Section  3,   hot-jever,  that  at  these  high  energies  the  i   =  0   approxima- 
tion gives  an  incorrect  result  for  the  inelastic  scattering  cross  section. 

1»22  Work  of  K.  Takayanagi 

The  next  advance  in  the  work  of  Zener  was  undertaken  by  K.  Takayanagi  ^    -• , 

The  higher  angular  momenta  are  taken  into  accoimt  by  the  approximation 

1 1/2 


k^-  Hi*  i)/r{|' 


where  R  is  some  siiitable  value  of  r  near  the  classical  turning  point  of  the  motion* 
One  then  replaces  k  in  the  radial  equation  by  Ic  and  has  only  an  s-wave  radial  equa- 
tion to  solve.  The  transition  probability  may  then  be  calciilated  from  the  result 
of  Jackson  and  Mott.  For  i    not  too  large,  the  gradient  of  J^  (j?  +  l)/r  is  small 
compared  with  the  gradient  of  the  potential  Ae"  ,  at  least  near  the  classical 
turning  point.  In  this  case  the  inclusion  of  the  slowly  varying  i ( {  +  l)/r 
term  in  the  wave  number  k  may  represent  some  reasonable  approximation  to  the  be- 
havior of  the  radial  equation  for  different  (small) j?  . 

K.  Takayanagi  and  S.  Kaneko '-  -^  have  applied  this  method  to  the  vibra- 
tional deactivation  of  0^  by  He.  ?or  large  i  ,  of  cotirse,  this  approximation  is 
unreasonable.  However,  the  authors  claim  that  the  major  contribution  to  the  cross 
section  comes  from  small  Jl    values.  In  the  transition  probability  calculated  by 
Jackson  and  Mott  the  main  factor  is  the  term  e~   f  ~  o  »  This  term  does  decrease 
with  increasing  j2  . 

1.23  Work  of  VJidom  and  S.H.  Bauer 

r2i 

..    Widora  and  S.H.  Bauer'--'  have  investigated  the  vibrational  deactivation  of 
COp  on  collision  vdth  H_0  using  the  semi-classical  method  developed  by  Zener '-'-I, 
The  method  is  extended  so  that  the  transition  probability  for  vibrational  deactiva- 
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tion  can  be  calculated  as  a  function  of  impact  parameter,  the  classical  analogue 
of  the  quantum  mechanical  angular  momentum.  A  spherically  symmetric  potential 
must  be  tised  in  order  to  solve  the  classical  equations  of  motion  for  the  inter- 
nuclear  distance  as  a  function  of  time.  This  solution  is  then  used  to  calculate 
the  transition  probability  as  a  function  of  impact  parameter  at  variovis  energies. 
The  transition  probability  decreases  vjith  increasing  impact  parsuneter,  as  expected. 
Also,  many  values  of  impact  parameter  give  a  significant  contribution  to  the  trans- 
iton  probability.  The  final  result  is  quite  sensitive  to  the  precise  values  of  the 
parameters  in  the  unperturbed  Lennard- Jones  12-6  potential  and  those  in  the  pertur- 
bation potential,  Widora  and  S,K.  Bauer  find  that  small  variations  in  these  para- 
meters result  in  a  large  variation  in  the  inelastic  scattering  cross  section. 
This  sensitivity  of  the  inelastic  cross  section  to  the  precise  values  of  the  po- 
tentials has  been  observed  in  the  present  work  (Section  3) • 

l,2h    Work  of  E,  Bauer 

Recently  E»  Bauer'- J  has  attempted  to  take  into  account  in  a  physically- 
reasonable  way  the  high  order  angular  momenta  which  contribute  significantly  to 
the  cross  section.  The  n-umber  of  angrilar  laoraenta  J^  that  contribute  significantly 
to  the  cross  section  is  given  approximately  by  kd  where  k  is  the  wave  number  and 
d  the  range  of  the  potential.  Since  k  is  relatively  large  for  a  molecule  because 
of  the  large  mass,  this  will  be  a  number  of  the  order  10-100,  Since  it  is  a 
hopeless  task  to  integrate  the  radial  equation  for  many  values  of  angular  momentum 
J^ ,  the  problem  is  reformulated  using  an  impact  parameter  b,  rather  than  the  angular 
momentum  J^,  The  problem  of  summing  over  x    is  thus  reduced  to  an  integration  over 
b.  The  wave  functions  are  obtained  by  the  WKB  method. 

The  method  of  impact  parameter  is  well  suited  to  problems  of  this  type. 
Because  of  the  large  mass  of  the  particles,  the  deBroglie  wavelength  is  much  smaller 
than  the  range  of  the  intermolecular  potential.  Consequently  one  can  speak  of  a 
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well-defined  orbit  for  the  collision.  This  means  the  WKB  approximation  for  the 
wave  functions  will  be  a  good  approximation.  On  the  other  hand,  it  can  be  shown 
that  a  completely  classical  treatment  is  inadeq-uate.  Assuming  that  the  coupling 

between  the  internal  coordinate  p  and  the  external  coordinate  r  is  induced  by  a 

t 
perturbation  potential  V  ,  one  can  calculate  roughly  the  average  momentum  trans- 
ferred in  a  collision.  Aen  it  turns  out  that  6p6x  <  n  for  most  values  of  the 
energy  that  are  of  interest.  Hence  a  quantum  mechanical  treatment  is  called  for. 
The  combination  of  the  WKB  method  for  the  vjave  functions  and  the  impact  parameter 
formulation  instead  of  partial  wave  analysis  is  suited  to  the  classical  and  quantum 
mechanical  features  of  the  problem* 

The  case  of  an  IL-  Hp  collision  is  treated  beacuse  there  exist  rather  de- 

ri3i 

tailed  calculations  of  the  intermolecular  potential  for  this  system^  -■.  Ultimately 
the  matrix  element  of  the  perturbation  potential  between  the  initial  and  final 
states  must  be  calculated.  Since  both  the  initial  and  final  states  are  highly  os- 
cillatory functions,  this  integration  is  not  easy  to  carry  out.  Approximations 
must  be  made  in  order  to  estimate  the  matrix  element.  An  effective  range  of  inter- 
action is  defined  which  includes  the  turning  points  of  the  orbit.  The  potential 
is  then  taken  to  be  zero  outside  the  limits  of  integration,  which    now  extend 
only  over  the  effective  range  of  interaction,  and  is  replaced  by  an  average  value 
within  these  limits.  The  precise  numerical  value  of  the  effective  range  is  un- 
clear, yet  the  integral  must  be  a  very  sensitive  function  of  it.  These  approxima- 
tions leave  the  numerical  values  of  the  calculation  in  doubt. 
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1.3  Experimental  work  on  interraolecular  potentials 
1.31  Gas  kinetic  and  other  experimental  data 

Much  work  has  been  done  in  recent  years  in  fitting  a  potential  to  a 
given  moleciilar  system.  Two  or  more  adjustable  parameters  are  varied  so  as  to 
obtain  the  best  fit  vdth  available  data.  Most  work  has  been  done  with  a  central 
potential,  since  there  are  not  enough  data  on  angular  variation  of  properties  of 
molecular  S3rstems« 

Two  potential  forms  have  been  exhaustively  investigated:  the  Lennard- 
Jones  6~12 


(7)  V  =  :Xr'^-  nx'^^ 


and  the  exponential-six 
(8)         V 
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where  r    is  the  value  of  r  for  which  V  has  a  spurious  maximum.  -Por  the  rare 
max 

gases,  a  is  of  order  12-15,  r  and  r  from  3-h  A,  and  e  from  .C026  to  .0086  ev. 
•=    '  'mo 

Data  on  variovis  equilibritim  properties  of  gases  which  involve  the  inter- 
molecular  potential  have  been  used  for  potential  fitting.  These  are:  the  trans- 
port coefficients  of  gases j  the  second  virial  coefficient,  which  appears  in  the 
equation  of  statej  and  crystal  data  -  the  equation  for  the  heat  of  sublimation  at 
0  °K,  and  the  equation  for  the  crystal  lattice  spacing. 

There  are  two  objections  to  the  use  of  these  data  in  attempting  to  fit 
(7)  or  (8)  to  any  molecular  system:  first,  these  data  refer  to  equilibrium  pro- 
perties of  gases  and  can  only  give  information  about  the  potential  for  relatively 
large  values  of  separation  of  the  moleciiLes.  In  the  study  of  inelastic  collisions, 
however,  we  are  interested  in  the  potential  for  very  small  values  of  nuclear  separa- 
tion. Second,  these  equilibriiim  properties  are  not  calculated  directly  from  the 
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intermolecular  potential,  but  rather  indirectly  by  averaging  over  a  large  aggregate 
of  molecules  in  interaction.  This  averaging  procedure  will  introduce  errors  of  its 
oim,  and  as  a  result  the  constants  in  the  potential  will  not  be  uniquely  determined. 

The  transport  coefficients  of  a  gas  are  given  as  functions  of  the  elastic 
cross  section  in  the  Chapman-Enskog  theory  of  non-ideal  gases.  The   relations  between 
the  transport  properties  and  the  cross  sections  are  derived  ;d.th  the  folloifJlng  as- 
siunptions:  (l)  only  elastic  binary  collisions  occur j  (2)  classical  mechanics  is 
validj  the  effects  of  quantum  mechanics  are  assumed  to  be  insignificant.  In 
addition,  a  spherically  symmetric  potential  is  postulated  in  order  to  carry  out 
the  integrations.  The  assumption  of  binary  collisions  limits  the  method  to  gases 
of  moderately  low  density.   Ihe  neglect  of  inelastic  scattering  makes  it  in^jossible 
to  calculate^ for  example,  the  coefficient  of  thermal  conductivity  for  a  gas  that 
can  transport  energy  via  its  internal  degrees  of  freedom,  e.g.,  Hp.   The  use  of 
classical  mechanics  restricts  the  theory  to  gas  temperatures  high  enough  that 
quantum  effects  are  indeed  negligible.  These  limitations  are  not  unduly  restric- 
tive for  most  molecules  of  interest.  Only  in  the  case  of  He  or  Hp  at  low  tempera- 
ture does  the  theory  fail. 

The  equation  of  state  of  a  non-ideal  gas  at  moderately  low  density  is 
given  by  the  equation 

(9)  pV  =  nkr(l  +  I  +  4r  +  ..«  ) 

where  B  and  C  are  functions  of  the  temperature  T  and  of  the  intermolecular  potential. 
If  classical  n^chanics  is  used  in  deriving  (9) ,  then  B  and  C  are  given  by 

(10)  B  =  I  J  f(r^2)  ^'^ 

where  f (r)  =  1  -  e"  ^^^^        and  U(r)  is  the  interatomic  potential. 


Hp  IS  anomalous,  anyway. 


-  12  - 
For  a  spherically  symmetric  potential, 
(11)  B  =  2n  fr^d  -  e-^(^^/^)  dr  . 


2 
Since  B/V  is  much  greater  than  C/V  at  normal  pressures,  and  the  equation  for  B 

is  much  simpler  than  that  for  C,  most  use  has  been  made  of  data  for  B,  the  second 

virial  coefficient. 

VarioTis  difficulties  arose  in  the  atter^Dt  to  make  a  reasonable  fit  of 

the  Lennard- Jones  potential  (?)  to  some  simple  molecxilar  systems.  T.  Kihara'-  -^ 

evaluated  the  third  virial  coefficient  for  the  rare  gases  using  the  Lennard-Jones 

potential  and  found  that  the  potential  bowl  -was  too  narrow.  This  is  to  be  ex- 

-12 
pected,  since  the  repulsive  r    term  is  not  accurate  over  a  large  range  of  dis- 
tances. The  correct  potential  behaves  near  the  origin  like  r~  rather  than  like 
r"  ■  .  It  has  been  found  impossible  to  fit  the  data  on  viscosity'-  -^  for  Hp,  COp, 
and  He  over  the  coir^lete  range  of  temperatures  for  which  data  is  available.  The 
parameter  £  in  (7)  must  vary  with  temperature  if  the  data  is  to  be  fit  for  moderate 
and  high  temperatures  simultaneously.  Yntema  and  Schneider  ^  J  encountered  similar 
problems  in  fitting  the  second  virial  coefficient  data  of  He,  Amdur,  Ross  and 
Mason'-  -'  found  that  tvro  different  sets  of  values  of  (r  ,e)  were  necessary  in  order 
to  fit  the  second  virial  data  of  COp  over  the  entire  range  0  <  T  <  600°C. 

These  difficulties  have  led  to  the  investigation  of  the  more  realistic 
exponential-six  potential  (8).  The  exponential  form  is  known  to  be  a  better  ap- 
proximation to  the  repoilsive  part  of  the  intermolecular  potential  than  the  form  Xr"  • 

Ifason  and  Rice  have  fitted  exponential-six  and  Lennard-Jones  12-6 
potentials  to  He  and  Hp  data.  L  -'  Here,  the  quantxun  corrections  to  the  second 
virial  coefficient  and  coefficient  of  Tdscosity  are  significant,  but  unfortunately 
can  only  be  estimated.  The  results  seem  to  indicate  a  better  agreement  with  the 
experimental  data  on  He  for  the  exponential-six  than  for  the  Lennard-Jones  potential. 
The  exponential-six  result  fits  the  data  on  viscosity  very  well  at  high  temperature, 
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whereas  the  Lennard-Jones  potential  gives  too  small  a  result  for  temperatures 
greater  than  about  800  K.  This  result  with  the  Lennard-Jones  potential  was 

ri5i 

also  obtained  by  Whalley  and  Schneider.  •-  -'  The  thermal  conducti-idty  does  not 
fit  particularly  well  with  the  exponential-six  potential  although  it  is  better 
than  the  Lennard-Jones.  Similar  results  are  obtained  for  Hp,  except  that  the 
thermal  conducti-vity  can  not  be  fitted  with  the  classical  theory  used  here.  Some 
typical  results  for  the  He  and  Hp  are  shovm  in  Table  I, 


He 

Hp 

Exp-slx 

L.- J. (12-6) 

Exp-six       L.-J. (12- 

a 

12.U 

•  •  • 

lij.O 

•  •  • 

r  ,a 

m'  0 

5.93 

5.1i3 

6.32 

6.22 

eA,°K 

9.16 

10»22 

37.3 

37.0 

Table  I 
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For  the  rare  gases  Ne,  A^  Kr,  and  Xe,  Ifeson  and  Rice  •-  -•  used   crystal 
data  to  find  r  and  e  as  a  function  of  a,  and  then  fitted  a  to  data  on  viscosity 
and  second  viidal  coefficient.  The  results  are  in  good  agreement  with  experiment. 
However,  for  these  rare  gases  a  Lennard-Jones  potential  gave  almost  the  same  results. 
It  was  not  found  possible  to  obtain  a  good  fit  to  the  data  on  viscosity  with  either 
of  the  two  potentials,  Thla  result  was  aliso  obtained  by  Vlhalley  and  Schneider.  ^^J 

In  an  attempt  to  resolve  the  difficulties  in  fitting  viscosity  data, 
Whalley  and  Schneider  reconsider  the  problem  of  fitting  the  exponential-six  and 
Lennard-Jones  potentials  to  Kr,  A,  and  Xe,  They  reject  the  use  of  crystal  data 
for  the  following  reason.  The  theoretical  formulas  for  the  dependence  of  the 
crystal  lattice  spacing  and  heat  of  sublimation  on  the  interraolecular  potential 
are  both  derived  on  the  ass\mption  of  the  additivity  of  the  potential.  This  means 
that  the  total  potential  energy  of  an  assemblage  of  molecules  is  taken  as  the  sum 
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of  the  potential  energies  between  isolated  pairs  of  molecules.  This  assumption  wag 
investigated  by  T.  Kihara  L  ^J  and  Jansen  and  Dawson  >-   J  who  found  it  led  to  incor- 
rect conclusions  about  crystal  structure.  VJTialley  and  Schneider  reject  the  use  of 
crystal  data  and  use  only  second  virial  coefficient  data  in  order  to  evaluate  the 
parameters  in  exponential-six,  12-6,  and  9-6   potentials. 

A  significant  fact  emerges  from  the  calculations  of  Whalley  and  Schneideri 
All  the  potentials  used  fit  the  virial  data  well.  Cn  the  other  hand,  for  the  three 
potentials  investigated,  there  is  a  large  spread  in  the  values  of  e/k,  and  a 
spread  of  nearly  10%  in  the  values  of  r  .  For  argon,  s/k  ranges  from  89,6  K  for 
a  9-6  potential  to  131.5  °K  for  an  exponential-six  potential,  r  varies  ftom  I4.O8 
to  3«78  A.  Similar  results  are  obtained  for  Kr  and  Xe. 

Whallejr  and  Schneider  conclude  that  the  coefficient  of  viscosity  is  re- 
latively insensitive  to  the  exact  details  of  the  interaction  and  can  be  moderately 
well  predicted  by  a  single  choice  of  e  and  r  .  The  reason  is  that  the  main  part 
of  the  temperature  variation  of  the  coefficient  of  viscosity  is  proportional  to 
the  square  root  of  the  temperature,  vihereas  the  part  depending  on  the  details  of 
the  intermolecular  potential  is  a  relatively  slowly  varying  function  over  a  large 
range  of  temperature. 

It  thus  appears  that  it  is  not  possible  to  obtain  a  unique  determination 
of  the  parameters  for  an  intermolecular  potential  by  fitting  to  transport  and 
second  virial  coefficient  data.  The  averaging  procedure  that  is  necessary  in 
order  to  go  from  the  properties  of  an  isolated  molecular  pair  to  the  properties 
of  a  gas  also  smooths  out  the  details  of  the  interaction  of  the  isolated  molecular 
pair. 


-15- 

1,32  Scattering  experiments 

It  is  also  possible  to  fit  an  intermolecular  potential  by  use  of  scat- 
tering data.  This  would  seem  to  be  a  more  satisfactory  procedxire  than  fitting 
to  macroscopic  data  such  as  transport  coefficients.  Scattering  cross  sections 
reflect  the  interaction  of  individual  molecular  systems  directly,  and  no  averaging 
over  a  large  assemblage  of  interacting  systems  need  be  made.  This  reduces  the 
chances  for  error  as  well  as  increases  the  possibility  of  obtaining  a  better  deter- 
mination of  the  parameters. 

Unfortunately,  little  molecular  scattering  data  is  available.  This  is 
due  to  the  difficiilty  in  carrying  out  molecular  scattering  experiments  at  the  low 
energies  of  interest.  Some  scattering  data  have  been  collected  for  alkali  atom- 
rare  gas  collisions  because  of  the  convenience  and  accuracy  of  the  hot  wire  method 
of  detection  that  can  be  -used.  The  results  are  useful  in  estimating  the  van  der 
Waal's  constant,  if  an  asymptotic  potential  of  the  form  -C/r  is  used  to  fit  the 
scattering  data. 

Recently  Amdur  and  Mason,  in  a  series  of  papers,  have  reported  the  cross 
section  for  scattering  of  the  rare  gases  A-A^-  -I,  Ne-Ne  ^  -',  Kr-Kr '-  -',  and 

r2^  r26i 

He-He'-  <J>L  J  ±-^   the  energy  range  from  500  to  2000  ev.  They  have  been  able  to 
fit  a  potential  of  the  form  V(r)  =  K/r  to  the  data  obtained,  using  classical  scat- 
tering theory.  The  potentials  so  obtained,  which  fit  the  data  very  well,  are  valid 
in  the  region  0,5  A  <  r  <  2A,  approximately.  Representative  values  of  K  and  s  are 
shown  in  Table  H, 


K 

s 

Range  of  validity 

A 

1.36  10"^ 

8.33 

2,18  A  <  r  <  2,98  A 

N 

5.10-^° 

9.99 

1.76  A  <  r  <  2,32  A 

Kr 

2.55  10-1° 

5.)42 

2.I42  A  <  r  <  3. Ill  A 

He 

7.55  lo"-'-^ 

5.9ii 

1.27  A  <  r  <  1.59  A 

Table  II.  Constants  K  and  s  for  the  intermolecular  potential  ?  =  K/r  derived 
from  scattering  data.  (V  in  ergs,  r  in  angstroms.) 
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1»33   The  potentisLi.  between  two  hydrogen  molecules 

The  interaction  of  two  hydrogen  molecules  has  been  calculated  anal;>-ti- 
cally  by  means  of  a  variational  method.  Because  of  the  way  the  calcialation  is 
performed,  the  interaction  can  be  written  as  the  svm   of  four  terras:  the  valence 
energy  due  to  the  overlap  of  thevave  functions^  the  eledzostatic  quadrupole- 
quadrupole  interaction  energy^  and  the  two  dispersion  energy  temis  varying  as  the 
inverse  sixth  and  eighthpowers  of  the  nuclear  separation. 

The  valence  energy  is  the  most  difficult  part  to  calculate.  Very  nany 
exchange  integrals  must  be  calculated  whether  a  perturbation  or  variational  treat- 
ment is  used.  DeBoer*-  -*  used  a  first  order  perturbation  (Heitler-London)  method, 

ri3i 

but  neglected  all  exchange  integrals.  Evett  and  Margenau*-  -'  did  a  similar  calcu- 
lation but  included  estimates  on  the  exchange  integrals.  The  result  for  the  valence 
energy  is  a  complicated  function  of  the  orientation  of  the  two  molecules. 

The  long  range  energy  of  interaction  is  much  easier  to  calculate.  An 
accurate  calculation  of  the  quadmpole-quadrupole  interaction  was  made  by  James 
and  Coolodge'-  -I  \ising  a  thirteen  terra  wave  function.  Margenau^  -I  obtained  a 
good  calculation  of  the  inverse  sixth  power  dispersion  energy  averaged  over  all 
orientations  of  the  two  molecules.  The  contribution  of  the  inverse  eighihpower 
dispersion  energy  to  the  total  interaction  is  very  small,  and  has  been  calculated 
as  a  function  of  nuclear  separation  only. 

The  result  of  all  these  calculations  gives  an  energy  of  interaction 
averaged  over  all  molecular  orientations  wluch  is  in  good  agreement  with  the 

r29i 

experimental  energy  calculated  fron  compressibility  data  •-  -".  The  potential  has 

a  zero  at  5.7  a  and  a  minimum  at  6«li  a  ,  Jiuch  effort  has  been  directed  at  ob- 
o  0 

taining  reliable  values  of  the  energy  for  molecular  separation  from  5  to  10  a^. 
This  is  the  region  of  interest  for  low  energy  collisions  which  are  important  in 
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calculating  macroscopic  properties  of  gases  (transport  coefficients,  and  the 

second  virial  coefficient).  However,  for  inelastic  collisions  it  is  necessary 

to  know  the  potential  in  the  region  2,5  to  5  a  . 

We  use  the  potential  of  Evett  and  Margenau  in  the  way  E,  Bauer  did. 

A  I'forse  or  exponential-six  form  can  be  fitted  to  the  calculated  potential  by 

matching  the  zero  and  ndniraum  of  the  potential,  and  the  slope  at  V  =  0.  In 

this  way  we  obtain  a  reasonable  approximation  to  the  potential  in  the  important 

region  from  2  to  5  a  • 
^  o 

2,  Inelastic  scattering  for  Hp-  Hp 
2.1  Formulation 

2»11  The  perturbation  scheme  and  the  potentials L J 

Consider  a  molecule  AB  colliding  with  another  molecule  (or  atom)  C.  So 
far  as  the  nuclear  motion  is  concerned,  A,  B,  and  C  are  regarded  as  point  particles 
moving  in  a  potential  field  which  is  determined  by  the  orbital  electrons.  This  is 
the  adiabatic  approximation  discussed  in  Section  1,  p.  1. 

We  use  the  coordinates  of  (l)  •  -^^or  a  collision  between  two  hydrogen 
molecules,  A  =  B  =  K,,   and  C  =  2A  =  2Mu.  Hence  the  reduced  masses  of  (2)  become 

(12)  M^  =  iy2  J       Mq  =  ^  • 

We  shall  consider  only  collinear  collisions  in  the  present  work.  The 
introduction  of  non-colD-inear  collisions  would  greatly  complicate  the  calculation 
by  bringing  in  more  independent  variables.  Since  we  do  not  know  the  intermolecular 
potential  very  well,  we  shall  use  a  simple  model,  i.e.,  a  spherically  symmetric 
potential.  The  potential  is  given  by  (3) : 

(13)  V(?,^  =  Vj(?^,^  -  V°(?,?^)  +V'(?;?^5  t-t^)       . 
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We  assume  that  V  is  very  much  smaller  than  V  .  Its  specific  form 
will  be  discussed  below. 

Since  we  are  concerned  with  a  transition  from  the  ground  to  the  first 
excited  vibrational  state,  the  assumption  of  a  harmonic  oscillator  potential  for 
V. (r  ,p)  will  be  good.   Thus  we  take 

(U4)  V°(?^,p)  =  ^  o?{-^.^/     . 

This  form  for  the  potential  V.  leads  to  a  simplification  of  the  perturbation 
potential  V. 

We  use  the  same  scheme  for  the  perturbation  potential  that  has  been 
described  by  Zener '- -^  *  ^ -'  and  Bauer  L-l,  The  perturbation  potential  arises  from 
the  expansion  of  V  (r,p)  in  a  Taylor's  series  about  p  =  p  : 


(15)     V°(?,^  =  V°(?,^  )  .  (p.  p  )  M    .  l(p.  p  )2  -e     ^_^ 


P  PQ  V  ^  /  P=Po 

The  first  term  on  the  right  is  the  unperturbed  potential  which  appears  in  (13). 
All  the  terms  beyond  the  first  constitute  the  perturbation  potential.  However, 
for  harmonic  oscillator  wave  functions, the  matrix  element  of  (p-p  )  between  the 
m-th  cind  (in  +  l)-st  state  vanishes  unless  n  =  1,  Hence  all  the  terms  in  (l5)  of 
higher  order  than  (p-  p  )  give  no  contribution  to  the  inelastic  scattering  cross 
section.  The  effective  perturbation  potential  is  then 

(16)  V'(?;p^,  p-p^)   =  (p-p^)  ^]  . 

Thus  the  effect  of  the  harmonic  oscillator  approximation  to  V.  is  to 
make  the  perturbation  potential  a  separable  function  of  p  and  r.  This  will  con- 
siderably simplify  the  perturbation  scheme  to  be  adoptedo 

We  are  concerned  with  a  transition  from  the  ground  to  the  first  excited 
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vibrational  state.  No  rotational  states  are  excited,  so  that  conservation  of 
energy  takes  the  form 

2     ,    x2.2 


e     • 


(17)  ii(co^-co^)      =     (k^-kpV/2M^ 
For  the  first   vibrational  state  of  Hp: 

(18)  Hoir.-  CO  )      =     20.10"^  e^/a.     =     0.^  ev       . 

10  '0 

A  Morse  potential  of  the  form 


(19)         \=    \^ 


g-2a(r-x)  _  g-a(r-x)| 


is  substituted  into  V  (r,p)  which  is  taken  as 

(20)  V°(?,5)   =  V^(?,,)  .  V„(?g,) 

where 

This  gives 

(22)  V°(?,^  =  2U^[cosh  ap  e-2^(^-^>-  cosh  f  e-^^^"^^]   . 

We  write  this  in  the  alteinate  form 

(23)  V°(?,?^)     .    U°[e-2-('-<i'.  e-^t"-^'] 

o  r  ~l2  1         cosh  ap 

TJ       =     2U^ [cosh  ap  /2     /cosh  aps     d     =    x  +  --ln  r: ^-wr     . 

e  o[_  "^o'  J  '  '^o'  a        cosh  ap  /2 

The  constants  are  chosen  to  fit  the  best  data  of  Evett  and  Jfeirgenau  L^  -' ; 

(2I4)  U°  -    U.18  lO''^  I-    ,      a    =    .916  a"^  J        d    =     5.6?  a      . 

e  a      '  o  0 

o 

V  (r,p  )  has  the  general  shape  shown  in  figure  2,  with  a  zero  at  r  =  d  =  5,6?  a  , 
e    0  o 

and  a  minimum  of  -  r  U°=  -o003  ev  at  r  =  d  +  In  2/a  =  6«i42  a  . 

4  6  '0 
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Intermolecular  Potential  V*  (r,  p  ) 


Figure   2 
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From  (16)  the  perturbation  potential  is 


„i,^jw  j^jk^     ,    ,    L  .  .     -2a(r-x)   .  ,  ap  -a(r-x) 
V  (r|p  jp-p  )  =  (p-p^)U^a  2  sinh  ap  e       -  sinh  -f-  ^ 


(25) 

which  can  be  written  as 

(26)      V'(?j  p^lP-iS-^)  -  (p.p^)u'a[e-2^('--'-  e-^'-^'] 

where 

'o'Ol 


u'a(p-p^)^T  =  5.86  10"^  ^/^'>  2  =  6.87  a   . 


Here  (p-p  )  -  is  the  matrix  element  of  (p~p  )  between  the  ground  and  first  ex- 
cited vibrational  states. 

2,12  The  wave  equation 

With  the  coordinates  (1)  described  above,  the  motion  of  the  center  of 
mass  separates  out  and  the  total  Hamiltonian  becomes 

(27)      H(?.?)   =  -^^l-t-  t*  ^l(?o.P>  *  0?'Pc>  *  ^'(?J?o'H>   • 

1  ^    e 

Since  the  perturbation  potential  V  is  sparable  in  p  and  r  as  we  see 
from  (26),  we  can  separate  the  Schrodinger  equation  by  writing  the  eigenstates  of 
the  unperturbed  sj-^tem  as: 

(28)  T   =2;   ^At)  Y.(^ 

where  0.  and  T.  are  solutions  of 
J     J 

¥e  treat  the  0's  as  known  harmonic  oscillator  wave  functions,  and  the 
potential  V  as  a  perturbation.  Our  problem  is  to  solve  (30)  for  the  T.(r)  and 
calculate  the  inelastic  cross  section  with  these  solutionso 
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Let  us  expand  the  functions  T.(r)  in  a  series  of  spherical  harmonics; 


(31) 


Y,(?)  = 


^ 


(2i+  DP,  (cos  ej  %,   (r)/r      j  =  0,  f 


'r  '1 


Substituting  into  (3C)  we  have 


(32) 


d^X 


L  . 


dr 


2M 


f :(?^Po^}  - 


r 


Xj   -   0 


We  use  perturbation  theory  in  order  to  obtain  the  inelastic  cross  section. 
Given  the  Hamiltonian 


and 


H  =  H  +  H    where  H  «  H 
o  o 


(H  -  E  )  u^(r)   =  0  , 
o   n   n 


first  order  perturbation  theory  gives  for  the  transition  probability  p_  that  the 
system  will  be  found  in  a  state  _f  with  eigenfunction  u„ 


2n 


Pfo=  T     Pfol  Pf 


•  i2 


where  p-dE^  is  the  number  of  states  in  a  range  dE^  about  E^,  and 


H, 


.^(r)   H  (r)   u^(r)  dt      . 


'fo    '    j     ^i 
Following  Schiffi^-'   we  obtain  for  the  cross  section 


2  k 


(33) 


d<r(a,p)     =       -^)     ^  )H^^p  d(cos  a)  dp 


,2riti' 


f    1    •    i2 

O 


In  order  to  determine  the  normaLization  of  the  functions  %^  (r),  we 
shall  calculate  the  cross  section  from  (33)  in  the  Bern  approximation,  in  which 
case  the  solution  T.(r)  of  (30)  is  just  a  plane  wave  for  both  the  initial  and 
final  states.  We  write  these  plane  waves  in  a  series  of  spherical  harmonics, 
using  the  well  knovm  expansion 
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(31)  e  °     =  T:  (2  5  +  l)i*  ^  (cos  e^)  j^  (k^r) 


(35) 


ik  r 


f  =  V  Co  {  +  ^^4 


T;;  (2^  l)i*Pj^(cos  e^jj^  k^)   . 


Here  we  take  the  z-axis  along  k  .  6  is  the  angle  between  k  and  r,  and 


©-  the  angle  between  k„  and  r.  Then 
1  1 


(36) 

vrhere 


Hf o  =    5^  5!^  ^2^  -"^^^2  J^  '+  l)r  "   P^  (cos  «^)P^,(cos  epd(cose^)djZ(j^^ 


(37)  V^^  =  j   jj^(k^r)v'(r)j^(k^)r2dr  . 

From  the  addition  theorem  for  spherical  harmonics, 

(38)  P„,(cos  e.)  =  21  %=-^  f?,(cos  e  )l^,(cos  a)  cos  ir^ 

where 

cos  a  =  cos  Q     cos  6-  +  sin  6  sin  ©_  cos  6 
0      1        0      1     '^ 

i.e.,  (a,p)  is  the  angle  between  k  and  k„« 

In  the  integration  over  0  in  (36)  only  m  =  0  from  (38)  will  contribute. 
But  then  the  integral  over  Q     has  the  term  •<P.  (cos  9  )Pg,(cos  Q   )>  S,(cos  a). 
This  integral  also  vanishes  unless  ^   =  \  }   then  (36)  becomes 


»fc 

(39) 

1 

H. 
fo 

From  (36) 

j  y~     (2£ +l)^P^(cos  ©^)P^(cos  a)d(cos  e^)^/].^ 


Un(2  H  +1)  P,  (cos  a)  v' 


dCT'(a,0)  =  j^  /-^j  H   Zi    l6n^(2X+l)(2]?'+l)P^(cos  a)P,(c<ns  a)|v^^|^ 
or     •■  N  / 

since  only  J^  =  X  contributes  because  of  the  orthogonality  of  the  P's« 
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The  normalization  here  is  determined  by  the  asymptotic  form  of  the 
spherical  Bessel  fvuictions,  jn  ,  namely 

(I4I)  i^   (kr)  ^^^     sin[kr  -  in/2]/kr 

and  hence  the  plane  waves  of  (3)j)  and  (3$)  have  the  asyn^jtotic  form 

(ii2)         e^^  -*  ZI  (2  le  +l)i^  P^  (cos  ©)  sin[kr  -  J^n/2]/kr  . 

The  Born  approximation  is  exact  only  if  V  -  0,  and  V  sufficiently  small. 
With  a  non-zero  potential  V  (r,p  )  present,  (31)  is  the  correct  solution  of  the 

6     O 

Schrodinger  equation.  Then  the  cross  section  derived  above  in  the  Bom  approxima- 
tion holds  also  for  the  case  V  /  C,  with 

(U3)        V^^  =  I  %^(k^r)  v'(r)  7^(k^)dr 

instead  of  (37) •  Now  the  potential  for  the  Hp-  Hp  system  goes  to  zero  for  large 
r.  Hence  asymptotically  (31)  must  approach  (U2).  Ihtis  the  normalization  of  the 
'p fl   's,  which  are  the  solutions  of  (32)^ is  determined  by 

where  the  phase  sliift  5.   represents  the  effect  of  the  potential. 

It  will  be  noted  that,(l40),  the  cross  section  in  the  Born  approximation, 
expressed  as  a  sum  over  phases,  reduces  to  the  \isual  form  of  the  Bom  approxima- 
tion if  we  carry  out  the  sums  and  integrations.  The  usual  form  of  the  cross 
section  in  the  Bom  approximation  is 

/°  2 


sin  Kr  „'/   ^  2, 
— ^ — V  (r)r  dr 


(h5)  d<r(a)     =   I^Y    ^ 

^22  ' 

vdiere  IT"  k  +  k.-  2k  k^cos  a.     From  (37)  we  have  for  V„ 

o      I         o  f  1 


d(cos  a) 


^fo  '  ^^  (V^^'^^^JA   (V)r2dr 
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Substituting  this  into  (39) 


(U6) 


r 


«fo=  ^" 


r  dr  Zl  (2|+l)j  (k^r)V  (r)j^  (V^^j^  (°°^  ''^  * 


The  sum  can  be  written  in  closed  form  since 


[31] 


sin  Kr 


T     (2J!+l)j^  (V)J^  (k^)P^  (cos  a)  -  ^^ 


Substituting  this  into  (li6),  we  have 


Hfo°     '' 


■( 


2^    irV   ^     sin  Kr 
r  dr  V  (r;     — ^^ — 


Than  the  cross  section  (33)  becomes 


d(r(a,p)     =     I41 


2^7 


J     r^dr  v' 


(r) 


sin  Kr 
Kr 


d(cos  a)dp 


Integrating  out  p. 


r 


2,    -n^'/   N     sin  Kr 
r  dr  V  (r)         ^ 


d(cos   a) 


Tjhich  is  just  (U5). 


Sujnmary 

We  now  have  the  following  scheme:  the  cross  section  for  inelastic  scat- 
tering is  given  by 


(I47) 

where 

(18) 

and 

(U9) 


(50) 


r  (r)  ^     i    sinrk.r  +  6,  -  in/2] 

^7  ^   '  r-»  00       k  |_  J  Jj?         'J 

The  functions    ^j^  are  the  solutions  of  the  equation 


J  =  0,f 


d^%, 


dr 


.2     ^e^o,         X       J((2 
V  -^  ^(^'Po>   -  — 


^} 


»     0     • 
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If  we  include  the  factors  l/k.  in  the  cross  section  rather  than  in 


the  Xg  i   we  have 


(51) 


(52) 


Im  (2\i 


We  convert  to  dimensionless  units  by  substituting  x  =  kr^  then 
(53)        --f-  + 


where 


dx 


^  Hi^+l)"' 


_    J 


2^2 


W"^  =  2M  U°  d  hi 


e  e 


:^^  =  0     j  =  o,f 


'i  -  ^j^ 


M°(x,p)=  e-2^'(^-Sj)_^-a.(x-Sj) 


e  ■"  o 


a«  =  a/k. 


With  these  dimensionless  units,  (53.)  can  be  integrated  on  the  Univac. 


2.2  Integration  of  the  radial  equation  on  the  Univac 

We  solved  (53)  by  a  modification  of  the  Runge-Kutta  integration  scheme 
by  S,  GiiiU'  J  which  is  particularly  suitable  for  high  speed  digital  computers. 
This  scheme  amounts  to  a  Taylor's  series  expension  about  a  point  x,  and.  in  effect 

calculates  the  function  and  its  derivative  at  a  point  x  +  h  from  the  function  and 

5  (v) 

derivative  at  x.  The  error  of  this  scheme  is  of  the  order  h  f ^   at  any  point, 

where  h  is  the  intei*val  of  integration  and  f    is  the  fifth  derivative.  Since 
max  Xq         =     1,  the  error  is  h  at  most.  A  s\ifficiently  small  interval  h  was 
tised  so  that  the  error  was  less  than  5  *  10  • 
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To  obtain  a  starting  value  for  the  integration,  we  •use  the  WKB  solution 


of  (53)  as  follovrs.     Let  x    be  the  solution  of 


[33] 


«^  mO  (^^^^,  .  Uj^ 


(5U)         •"■  "  '7  "e  ^^'Po^  "  ? 

0 


=  0 


For  J^  »  0,  X  is  equivalent  to  the  classical  turning  point  of  the 
motion  in  the  potential  given  hy   (23) •  Then  x  <  x  j  the  VJKB  solution  of  (53)  is 


(55) 

where 


^WKB 


|q-V2! 


exp 


Q  dx 


X  <  X 


W 


The  solution  in  the  region  x  >  x  is  given  by 


(57) 


Q~  '  cos 


Q  dx  - 


n 


X  >  X   • 

0 


The  WKB  solution  (55)   was  calcvilated  beginning  from  x  =  x     and  proceeding 
backwards  towards  x  =  0  (i.e.,  r  =  O).     When  the  solution  so  obtained  is  less  than 
10     ,  we  have  reached  the  asymptotic  region  in  which  the  WKB  solution  is  valid. 
The  starting  values  for  the  integration  are    ^  „    and  /J:^    j,    at  the  first  value  of 
X,  X  =  X. ,  at  which   /C^ja,  ^10     .     This  assures  us  of  not  having  neglected  any- 
significant  contributions  to  V  „     in  (1^8)   for  x  <  x.  • 

We  now  integrate  {^3)   forward  using  the  Runge-Kutta-Gill  scheme.     At 
the  turning  point,  x  =  x  ,  we  compare  the  ratio  of    "X  /%    (■^^ere    %  is  the  solution 
of  (53))  with  that  predicted  by  the  WKB-Langer  approximation;'--^  -' 


(58) 


C  - 


L^ 


dx/ 
/x=x 
/   c 


^^yiAb        r(i/3) 

The  two  ratios  disagree  by  about  2-3^o,  but  this  can  probably  be  attributed  to 
the  inaccuracies  of  the  approximation  used  to  obtain  (58),  rather  than  to  any  in- 
accuracy in  the  Univac  integration. 
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In  the  region  r  -»oo,  we  know  that  the  asymptotic  form  of  the  solution 
of  (53)  is  given  by  (52).  The  asymptotic  region  is  reached  when  the  roots  of  %j^ 
occur  at  a  distance  n  apart.  When  this  happens,  the  Univac  program  stops  integra- 
ting the  equation  (53),  which  is  relatively  time-consuming,  and  produces  the 
simple  sine  wave  (5l)  with  a  proportionality  constant  determined  by  the  preceding 
values  of  % ^    at  its  peaks.  The  integration  is  ended  at  some  fixed  value  of  r 
(12  a  )  for  ^ich  the  value  of  V  in  (26)  is  siofficiently  small  so  that  there  will 
be  no  significant  contribution  to  V   in  (ii8)  from  lar gar"  values  of  r. 

As  a  check  on  the  forward  integration,  the  WKB  solution,  (57)  was  calcu- 
lated simtiltaneously.  After  the  asymptotic  region  of  the  WKB  solution  is  reached, 
the  values  of  x  at  the  nodes  and  peaks  of  the  WKB  and  R-K-G  solutions  differed 
by  less  than  10  .  For  a  second  check,  the  solution  Xfl  w^s  substituted  back 
into  the  differentia],  equation  (53).  A  simple  seven-point  formula  for  the  second 
derivative  was  used.  The  solution  %  ^    was  found  to  satisfy  the  differential  eqxia- 
tion  to  i-rithin  an  error  of  the  order  of  3  to  5*10   at  any  point. 

The  calculation  of  the  matrix  elements  V   in  (U8)  was  carried  out  in  a 
perfectly  straightforward  manner.  The  interval  of  integration  h  in  the  RKG  inte- 
gration scheme  for  the  radial  equation  was  so  chosen  that  both  %»    ^^^^   ^^ 
"Xo  (kJ')  in  (ii8)  were  known  for  exactly  the  same  values  of  r.  This  avoided  any 
interpolation  in  the  radial  functions  %^    .  Hence  V^^  is  calculated  to  exactly 
the  same  accxiracy  to  which  the  %n     are  known. 

The  total  cross  section  is  then  found  by  integrating  (51).  A  comparison 
with  a  summation  rather  than  an  integration  over  X    showed  that  the  partial  cross 
sections  C.     =  (2^+1)  W^^     are  a  sufficiently  smooth  function  of  i    so  that 
the  integration  was  Just  as  acciirate  as  the  summation. 
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2,3  Resiilts 

2.31  Results  of  Una vac  integration 

The  potentials  used  and  typical  radial  functions  are  shown  in  Figure  3» 
The  corrected  result  for  the  cross  section  obtained  on  the  Univac  is  shown  in 
Figure  U»     The  results  obtained  by  E,  Bauer  are  also  shown.  The  correction  to 
the  Univac  cross  section  comes  about  as  follows. 

The  original  results  obtained  using  (5l)  are  not  flxix-conserving  at 
energies  greater  than  about  5  ev.  That  is,  the  partial  cross  section  of  ($1) 


(^''  -*  =  rb  [Br  '''*''  Kof 


[35l 

exceeds  the  maxiinum  allowable  value  of  the  partial  cross  section*^  ••, 

(60)  max  cf     =    ~~    ^^  ^  *^^      ' 

^  of 

The  partial  cross  section  of  (60)   is  tlte  largest  possible  if  the  flux  of  particles 

ig  to  be  conserved. 

The  violation  of  flux  conservation  at  these  high  energies  comes  about 

because  the  molecules  approach  very  closely  at  such  energies.     Then  the  wave 

t    o 

functions  will  be  appreciably  large  in  the  region  where  V  >  V  ,  i.e.,  near  the 

origin.   But  the  -perturbation  scheme  iised  here  is  a  reasonable  approximation 

only  as  long  as  V  «  V^. 

e 

From  a  formal  point  of  view,  the  violation  of  flux  conservation  arises 
because  the  scattering  matrix  of  the  perturbation  scheme  is  not  unitaiy.  This 
can  be  corrected  by  working  with  the  reaction  matrix  K  instead  of  T  =  S  -  1.  The 
scattering  matrix  S  is  given  by  1  -  ^  iK/1  +  ^  iK,  so  that  S  is  unitary  so  long  as 
K  is  hermitian  ^  J.  we  have  taken  V   to  correspond  to  the  eigenvalties  of  the 

reaction  matrix  K,  Since  these  eigenvalues  behave  as  tan  £«  ,  we  replace  Iv  I 

I  '  |2/     I  '  i2  2  2 

t>y  V^^  /  -^  *    rfo  '  i«®»>  replace  tan  5.  by  sin  5*  • 
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E(ev) 


2.5    3.0  4.0  50  6.0 

Distance  r  (Oo) 

Potentials  and  Initial  and  final  wave   functions 
for    Eo  =  l.l  ev,  i=0 


Figure  3 
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Then 
(61) 


Ivkl' 


*         of  1  +  v^  r 


fol 

It  is  this  corrected  cross  section  which  is  shown  in  Figure  L(.  At 
the  energies  where  this  correction  is  needed,  however,  the  inelastic  cross  section 
calciJ.ated  in  the  present  work  is  entirely  incorrect.  At  those  energies  the  pro- 
bability of  electronic  transitions,  i^iich  we  have  neglected,  becomes  significant. 
The  inelastic  cross  section  calculated  here  is  consequently  physically  unrealistic 
for  energies  greater  than  5  ev. 

It  is  clear  from  Figure  h   that  there  is  disagreement  vjith  the  results 
of  E.  Bauer  so  far  as  the  energy  dependence  of  the  cross  section  ia  concerned. 
The  reasons  for  this  have  already  been  pointed  out  in  Section  1,  ^  I.2I4,  But 
we  do  verify  here  that  one  has  to  sum  over  many  partial  waves.  This  is  shoim  in 
Figure  ^  where  the  partial  cross  section  is  plotted  as  a  function  of  a  for  both 
the  present  calculation  and  for  the  results  of  E.  Bauer.  In  Fi,gure  6  we  have 
plotted  on  a  linear  scale  the  cross  section  obtained  by  the  Univac  and  that  of 
E.  Bauer  in  the  energy  range  near iiires hold,  0»5  to  2  ev. 

To  obtain  checks  on  the  energy  dependence  of  the  inelastic  cross  section, 
and  in  particular  on  the  expected  value  of  energy  for  which  the  cross  section  has 
a  maximuiii,  various  approxLinate  solutions  to  the  problem  have  been  considered. 
These  are  discussed  in  the  following  section.  AH  these  approximations  indicate 
that  the  energy  dependence  of  the  inelastic  cross  section  in  Figure  U   is  physi- 
cally reasonable  and  that  the  maximum  can  be  expected  to  occur  at  a  relatively 
high  energy. 


53 


fo 

Different, 

arbitrary 

units 


Partial  cross-section  en  (X)  as  a  function   of  angular  momentum  Z 

for  various  energies 

Figure   5 
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Eo  (ev) 

Total   cross- section!  comparison  of  UNIVAC  results  with 
Bauer  (1955)  for  energies  near  the  threshold 


Figure  6 
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2.32   Approximations 

2.321  Resonance  approximation 

The  intermolecular  potential  can  be  looked  upon  as  a  dist-urbing  force 
applied  to  a  simple  harmonic  oscillator  of  frequenc;;-  v,  which  is  the  molecule 
that  can  be  excited  to  the  first  vibrational  state.  In  order  for  any  appreciable 
forced  oscillation  to  be  produced,  the  time  collision,  f  ,  should  not  be  too 
great  compared  with  the  natural  period  l/v.  Hence 

(62)  f  V  »  1      weak  excitation  , 

fv  '^  1      maxiraara  excitation  • 

Since  t^d/v  where  d  is  the  range  of  the  intermolecular  potential  and 
V  is  the  relative  velocity  of  the  colliding  molecules,  the  condition  for  raaximura 
excitation  becomes 

(63)  dv/v  -'  1  . 

Writing  v  =  AE/Ii,  we  obtain  k  '^30  a"       for  a  maxiiuum  cross  section, 
or  E  ~'   7  ev.      This  energy  is  of  the  same  order  of  magnitude  as  the  actual  maximum 
energy,   25  ev» 

2.322  Bom  approximation 

We  calculate  the  cross  section  for  inelastic  scattering  in  the  Born  ap- 
proximation vising  (U5)  for  the  cross  section  and  (26)  for  the  perturbation  potential. 
The  integrations  can  be  carried  out  simply.  The  results  are  sho^m  in  Sigure  ?• 
It  is  seen  that  the  maximum  in  the  curve  occurs  at  an  initial  energy  of  about  7  ev. 
As  an  independent  check  on  this  result,  the  cross  section  can  be  calculated  from 
the  alternate  form  of  the  Bom  approximation,  (UO),  as  a  sum  over  phases.  V   is 
given  by  (37)  and  involves  for  the  case  of  the  Morse  potential  the  integral 
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00 


where 


therefore 


I  =  r   e-^  j^  (k^r)j^  (k^)r^dr 


i^   (p)  =  ("/2p)^/2  j^^_^^^  (p) 


I  = 


,00 

n      f    -ar 


2(k^kp-^/^  ;o 


i;^       ^     Jji-.iA^V^'^^+iA^V^^  ^- 


of  the  second  kind 


This  can  be  written  in  terms  of  Qp  ,  the  associated Legendre  function 
.-..[37] 


2(k^kp 

where  the  prime  denotes  differentiation  -with  respect  to  the  argument.  An  ex- 
pression for  the  cross  section  is  obtained  in  terms  of  a  sum  over  H  of  these 
Legendre  functions.  This  sura  was  evaluated  for  two  energies,  and  the  answers 
agreed  with  the  simpler  form  cf  the  Born  approximation. 

The  Bom  approximation  result  in  Figure  7  is  not  flux-conserving  over 
the  entire  energy  range  of  interest,  that  is,  up  to  an  energy  so  large  that 
electronic  transitions  become  significant.  Furtliermore,  it  is  larger  by  several 
orders  of  magnitude  at  every  value  of  energy  than  the  cross  section  of  Figure  U» 
Consequently,  the  Born  approximation  is  of  absolutely  no  use  in  calculating  an 
inelastic  cross  section.  However,  as  a  mathematical  approximation  it  does  indi- 
cate the  expected  value  of  the  energy  for  the  maximum  cross  section. 

2.323  High  energy  approximation 

The  radicil  equation  we  must  solve  is  given  by  (53).  It  will  be  noted 
that  if  we  approximate  V  (r)  by  a/r  ,  where  a  is  a  positive  constant,  the  solutions 

6 

Xg    become  spherical  Bessel  functions  of  non-simple  order.  Furthermore,  if  we 
also  approximate  V  (r)  by  p/r  then  V„  in  (I48)  can  be  written  in  closed  form. 
In  short,  the  cross  section  can  be  calculated  analytically  in  this  approximation. 
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This  replacement  of  V  (r)  by  a/r  -will  surely  not  be  good  except  at  very  high 
energies  with  initial  energy  greater  than  about  2  ev.  However,  this  will  give 
us  an  indication  of  the  high  energy  behavior  of  the  cross  section. 
Substituting  V°(r)  =  a/r^  into  (53) 


e 


dx      I       X 


%t        =  0 


O  0  0 

where  a'  =  2ma/  'h  .   (a  must  have  dimensions  energy-cm  in  order  for  a/r  to  have 
dimension  energy,  as  it  should.  Hence  a'  is  dimensionless.)  The  solution  to  (6I4) 
is 

(65)  ;}^j^  =  j^^i/2  ^""^ 

which  has  the  correct  asymptotic- form  (ii9)  •  Here 
(66)         v(v+l)  -   S(jf +1)  +  a' 

2v  +  1  =  2[({  +  1/2)^+  a']^/^  . 

^^^^  ^w^a.v,^.v.  ^v^wv  ^x^  vv^v^/  ,,^,^<,  ^v.   ^.^.^v-xx      J  in  order  that    Xn     be  finite  at  r  =  0» 
Substituting  this  solution  into   (U8),  we  have 


r  dr 


2(k^kp 


^     /^    ?  Vi/2  (V>  Vi/2  (V) 


This  integral  can  be  written  in  the  closed  form 

Hence 

,v+l/2 


4  =  j^  ^  (''A'' 


(68)         (T  = 
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Then  the  cross  section  in  (U?)  is  just 

2  hn     rr-     /o()  .-,>  f,.    hr   ^2v+l 


^  '    o 


On  comparing  the  values  of  the  Morse  potential  of  (23)  with  the  approxi- 
mation  l/r  ,  it  is  seen  that  a  value  of  a'  of  about  5000  will  give  a  fairly  good 
correspondence  between  a/r  and  V  (r).  This  value  of  a'  will  be  good  at  least 
for  energies  high  enough  that  the  turning  point  is  not  greater  than  about  2,2a  , 
This  will  be  true  for  energies  greater  than  about  3  ev.  The  inelastic  cross 
section  was  calculated  from  (68)  with  three  different  values  of  a';  5000,  3000, 
and  7000.  These  are  shown  in  Figure  8.  It  is  clear  that  the  large  variation  in 
the  values  of  a'  produces  only  a  slight  change  in  the  position  of  the  peak  of  the 
curve.  The  value  of  the  cross  section  at  the  peak,  hovrever,  varies  greatly  with 
the  change  in  a'. 

As  a  check  on  the  accuracy  of  the  cross  section  in  this  apprcKiraation, 
we  can  compare  (68)  with  the  cross  section  obtained  from  the  Born  approximation 

for  a  p/r  potential.  In  this  case  the  cross  section  is  given  by  {\x\) •     Since 

to  f 

we  are  using  V  =  p/r  ,   the  integration  for  V       is  carried  out  like   (6?)   with  the 

only  difference  that  we  have     ^   now  instead  of  v: 

|2     In    ^       1  „,   /^^2     +1 

o 

We  see  from  (66)  for  (2v+l)  that  if  j^  is  large  enough  compared  with  a' 
that  we  can  neglect  a',  then 

(2V+1)  =  2[(j^+l/2)%  a'J    -j^rt,   2(1^+1/2)   -  2  J2 +1  . 
In  this  case  (68)  reduces  to  ^^9)  y   as  it  should  in  the  limit  of  high  energies. 
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2.33  Low  energy  cross  section 

The  low  energy  dependence  of  the  inelastic  cross  section  on  energy  is 
of  some  interest  in  the  present  calculation.  It  is  in  this  region  that  the  s-ijave 
approximation  may  give  reasonable  results.  Although  it  is  known  that  this  approxi- 
mation will  not  do  for  a  problem  in  inelastic  scattering  because  of  the  importance 
of  the  higher  phases,  many  theoretical  calculations  have  been  done  for  the  case 
X    =  0   since  it  is  possible  to  handle  it  analytically  for  some  simple  forms  of 
the  interraolecular  potential.  We  shall  compare  one  calculation  at  a  very  low 
energy  (initial  energy  at  the  threshold,  0. ^U  ev)  vdth  previous  theoretical  calcu- 
lations in  the  s-wave  approximation, 

Zener  ^  -^  considered  the  potential 

(70)  V  =  D[e-2"^^-^>-2e-^(^-'°>]   . 

This  potential  has  a  minimum  of  -D  at  r  =  b.  It  is  a  slight  modification  of  the 
Morse  potential  used  in  the  present  work.  Zener  uses  semi-classical  methods  to 
find  the  transition  probability  for  a  head-on  collision  with  H.   =  0   since  there 
are  not  available  analytic  continuous  spectrum  solutions  to  the  radial  equation 
for  the  potential  (70).  The  result  is 


8„1;  ^(^^  Vi  1   1    C0sh^(q^-qJ(Vn) 

"^    a%     Mg(M^+  M^+  M^)'^       ^  -  ^     sinh''  n(q^-  q^) 


where 


(72)     cos  -^   =  (1  +  ^  )     e''  =  ^  m""  #    "^  =  I  ( V  \^ 

q^  =  m'  v^Aa  j  q^  =  rfv^  5   M^"  =  M^(>^-  M^)/M^+  F^+  M^  . 

In  the  present  work 

^A  =  '^^     \-     ^  =-  \>     %'     V'^  '  Sn  =  k^a  ,  M^  =  r^  =  M^ 
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Also,  n  =  0,  and  we  must  take  the  positive  sign  in  (71) •  Then 

It  is  easily  shown  that  this  formula  gives  an  incorrect  value  for  the 
inelastic  cross  section  when  the  energy  is  very  large,  VJhen  k  -»oo,  k  — >  k  =  k, 

Tf  0  0 

and  7  "^  7  •  ^Iso,  since  k  -  72  =  k„,  where  72  represents  the  energy  difference 

2 
2K,Zi^E/^  between  the  ground  and  first  excited  vibrational  states, 

(VV^W  =  72 


Hence 


where 


k^-  k^  ~»  72/2k  . 


2        2 
p„  — >  const,  cosh  A/k/sinh  2A/k 


A     =     le  n/a     , 

2  2  2 

For  k  very  large,   cosh    A/k  •— »   cosh     0  =  I5   sinh     2A/k  approaches 


(2AA)^ 

'Ihen 

,2 

(7I1) 

_^       k 

for  k  very  large. 

Since  the  partial  cross  section  is  just  o^     =  (2j^  +l)n/k  k    p     ,    (fU) 

gives  a  i«sult  for  the  cross  section  vrliich  is  independent  of  energy.     It  is  known, 

-1/2 

however,  that  the  cross  section  must  decrease  as  (E)  '  for  large  energy. 

At  very  low  energies  near  the  threshold,  however,  we  may  expect  (73)  to 
be  a  good  approximation.  We  have  calculated  the  inelastic  scattering  cross  section 
on  the  Univac  at  the  threshold  energy,  E  =  ,5)4  ev,  and  with  a  final  energy 
E^  =  7»10   evo  Ihe  same  cross  section  was  computed  from  (73) •  The  results  are 
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shown  below. 

Univac 1.26  •  lO"^^ 

Zener   1.31  *  lo"^^ 

Even  at  this  very  low  energy  the  Jc  =  0  approximation  gives  a  resvilt  that  is  too 
large  by  a  factor  of  10. 

3.  Sensitivity  of  the  cross  section  to  variations  of  the  potentials 

We  novT  examine  the  dependence  of  the  inelastic  cross  section  on  the 
parameters  of  the  unperturbed  and  perturbation  potentials.  Since  the  intermolecular 
potential  is  not  known  too  precisely  for  most  systems  of  interest,  it  is  important 
to  know  Just  how  the  inelastic  cross  section  varies  when  the  potentials  undergo 
small  changes.  The  potential  of  the  H^-  PL  system  u^ed  in  Section  2  has  three 
adjustable  parameters  f cr  which  values  have  been  chosen  to  fit  the  best  available 
data  for  the  H^-  Hp  interaction.  Since  these  data  are  uncertain,  especially  in 
the  important  region  where  the  potential  is  steeply  repulsive,  it  represents  at 
best  only  an  approximation  to  the  actual  potential  between  two  H^  molecules.  We 
shall  vary  two  of  the  parameters  in  this  potential  in  a  systematic  way  to  determine 
Just  how  sensitive  the  inelastic  cross  section  is  to  small  changes  of  these  parameters. 

There  are  two  distinct  types  of  variation  to  be  considered.  Fii'st,  with 
a  fixed  unperturbed  potential  V   (r,p  )  we  can  vary  the  parameters  in  the  perturba- 
tion potential.  Since  the  perturbation  potential  is  derived  in  an  approximate 
fashion  from  the  unperturbed  potential,  a  knowledge  of  the  variation  of  the  in- 
elastic cross  section  viith  small  changes  in  the  perturbation  potential  is  necessary. 
From  it  we  can  see  how  accurately  the  perturbation  potential  must  be  known  in  order 
to  obtain  a  given  accuracy  in  the  inelastic  cross  section. 
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Second,  we  can  vary  the  parameters  in  the  -unpertiirbed  potential  V  (r,p  ). 
We  are  particularly  interested  here  in  the  sensiti%d.ty  of  the  cross  section  to 
variations  in  the  steepness  of  the  repulsive  portion  of  the  potential, 

3»1  Sensitivity  to  the  perturbation  potential 
The  perturbation  potential  is  given  by 

(75)  M?,-;„iP-P„)  -  (?-?„)  \^J^__^       . 

From  Section  2,  (26)  and  (29),  the  potentials  are 

(76)  V°(?,p^)  =  U  Je-2^(^-^^-  e-^^^-^^] 

(77)  V'(?,^^,p-^^)  =  U'a(p-p^)  [e-2a(r-z)  _  ^-a(r.z)J 

U^  =     I4.I8  »  10"^  '  eVa^j     u'a(p-p^)Q^  =     $.86'10"^  ®^/%J  ^  =  '^^^  ^'^ 

z     =     6.87  a  d     =     5»67  a 

o  0 

Ihe  wave  functions  calculated  in  Section  2  are  tised  here,  but  the  para- 
meters a  and  d  are  now  varied  about  the  'base'  values  of  (77) •  The  cross  section 
as  a  function  of  energj--  is  calculated  for  each  new  value  of  a  and  d. 

The  results  of  varying  the  steepness  parameter  a  are  shown  in  Figure  9« 
Two  values  of  a  roughly  5/o  greater  and  smaller  than  the  base  value  of  ♦916  a' 
were  chosen.  We  have  also  shown  the  result  previously  obtained  in  Section  2  using 
the  potentials  (76)  and  (77) •  We  see  that  the  cross  section  is  a  very  sensitive 
function  of  the  parameter  a,  A  change  of  roughly  y'/o   causes  a  variation  in  the 
cross  section  by  a  factor  of  about  2«5,  This  marked  sensitivity  to  the  value  of 
the  parameter  a  is  to  be  expected  from  the  following  considerations. 

Referring  to  Figure  3  in  Section  2  we  see  that  the  turning  point  ^f  the 
classical  motion  for  the  initial  and  final  states  at  typical  values  of  the  energy 
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occitr  in  a  region  where  the  perturbation  potential  is  a  very  sharply  varying 
fimction  of  r.  The  inelastic  cross  section  is  a  function  of  the  integral 

I  T^(r)v'(r)Yj^(r)r^dr  . 

In  the  region  near  the  turning  points  of  the  classical  motion  the  product 
Y  (r)  Y^(r)  is  a  quickly  varying  function  of  distance.  Now  the  most  important 
contribution  to  the  integral  comes  from  the  first  peak  in  the  integrand  and  this 
first  peak  occurs  near  the  turning  points.  (This  can  be  seen  approximately  from 
the  MKB  representation  for  Y  (r)  and  Y„(r).)  Since  the  perturbation  potential  is 
also  a  steep  function  of  r  near  this  first  peak,  it  is  clear  why  the  inelastic 
cross  section  depends  critically  on  the  value  of  the  perturbation  potential  in 
the  region  0  <  r  <  r  (where  r  is  the  value  for  which  V  is  zero) .  The  perturba- 
tion potential  in  this  region  is  in  turn  a  very  sensitive  function  of  the  parameter 
a,  as  can  be  seen  from  (77) •  Consequently  the  variation  of  the  cross  section  with 
changes  in  a  is  very  marked. 

From  the  form  of  the  perturbation  potential,  (26),  it  is  clear  that  a 
change  in  the  parameter  d  will  produce  essentially  the  same  effect  as  changing  the 
parameter  a.  Changing  the  parameter  d  moves  the  whole  perturbation  potential 
bodily  along  the  r-axis.  In  the  region  0  <  r  <  r  ,  which  is  the  significant  re- 
gion in  this  problem,  this  is  equivalent  to  multiplying  the  potential  by  a  constant. 
However,  the  effect  on  the  perturbation  potential  of  changing  the  parameter  a  is 
the  same,  at  least  in  the  region  of  interest.  Hence  we  may  expect  quite  similar 
results  for  a  change  in  d  as  those  Just  described  for  changes  in  a.  For  this 
reason,  we  did  not  carry  out  the  calculations  for  variation  of  the  parameter  d» 
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3«2  Sensitivity  to  the  -unperturbed  potential 

We  consider  now  the  variation  of  the  inelastic  cross  section  with  changes 
in  the  parameters  of  the  unpertiirbed  potential.  It  is  evident  that  variations  of 
the  unperturbed  potential  will  produce  a  variation  in  phase  shift  of  the  radial 
functions.  Ihese  variations  in  the  phase  shift  of  the  initial  and  final  radial 
functions  will  cause  changes  in  the  cross  section  for  the  same  reasons  discussed 
in  S  3.1.  The  first  peak  in  the  product  T  (r)Y„(r)  will  be  displaced  along  the 
r-cixis  because  of  the  phase  shift.  Since  the  perturbation  potential  is  very 
quickly  varying  in  this  region,  the  cross  section  iiill  vary  sharply  with  the  varia- 
tion in  the  phase  shift. 

We  proceed  to  examine  this  effect  in  two  stages.  First,  we  calculate 
the  variation  of  the  cross  section  with  phase  shift  in  the  radial  functions.  That 
is,  we  deliberately  insert  a  phase  factor  5  in  the  expression  for  the  transition 
matrLx  element  V.  (Section  2,  (5l)): 

ID 


(78)         yI^   =   ]  T^(r  +  5)  V*(r)  f^(r)r^dr 


fo 


and  calculate  the  cross  section  as  a  function  of  6,  Second,  we  vary  the  parameters 
a  and  d  in  (76)  and  observe  the  change  in  5  and  consequently  in  the  cross  section. 

A  first  attempt  at  this  program  was  made  with  the  exp-six  potential 
(Section  1,  (8)): 

(79)         V  .  -^     I  2  axp|a(l.il)|  -  & 


a 

This  potential  form  is  of  interest  becavise  it  has  been  widely  used  in  recent  years 
with  moderate  success  for  calculating  transport  coefficients  and  other  equilibrium 
properties  of  gases.  It  has  a  spurious  maximum  at  a  small  value  of  r  of  the  order 
2  -  2.5  a  .  This  is  of  no  conseq-aence  in  calculating  the  equilibrium  properties 

0 

of  gases  since  they  depend  on  the  behavior  of  the  potential  near  V  =  0© 
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Values  of  e,  a  and  r  were  chosen  so  that  (79)  fits  the  available  data 
on  H^-  H^*  The  radial  functions  calculated  with  this  potential  agreed  very  closely 
with  those  calculated  from  the  Morse  potential,  (77) •  For  this  Morse  potential 
the  transition  probability  is  a  monotonic  decreasing  function  of  j?  .  The  transi- 
tion probability  for  the  exp-six  potential  (79),  however,  showed  a  wrong  functional 
dependence  on  the  angular  momentum  H  ,     It  increased  from  some  value  at  J^  =  0  to 
a  maxiitiun,  and  then  decreased.  Since  the  radial  functions  and  unpertiirbed  poten- 
tials were  essentially  the  same  for  both  the  Morse  and  exp-six  potentials,  we 
could  expect  the  transition  probability  to  be  the  same  function  of  H    for  both 
potentials • 

This  behavior  is  easily  accounted  for  by  the  spxirious  maximum  in  (79). 

Since  the  pertxirbation  potential  is  derived  from  the  unperturbed  potential  it  too 

exhibits  a  maximum  at  r  ~'  2  -  2,5  a  *  This,  however,  is  a  serious  defect  in  the 

o 

perturbation  potential  because  of  the  relatively  high  energies  we  must  consider. 
At  these  high  energies  the  turning  points  of  the  classical  motion  for  small  angular 
momentim  are  near  this  spurious  maximum  in  the  perturbation  potential.  Consequently 
the  first  and  most  important  peak  in  the  integrand  of  the  transition  probability 
occurs  in  a  region  where  the  perturbation  potential  is  too  small.  At  higher  values 
of  angular  momentum  the  turning  points  move  out  to  greater  values  of  r  where  the 
perturbation  potential  is  larger, and  behave  correctly  as  a  function  of  r.  This 
accoiints  for  the  observed  (incorrect)  maximum  in  the  dependence  of  the  transition 
probability  on  i  • 

For  this  reason  the  exp-six  potential  is  of  no  tise  for  calculating  in- 
elastic cross  sections  where  the  value  of  the  perturbation  potential  at  small  r  is 
important.  Although  this  latter  behavior  is  of  no  consequence  in  the  calculation 
of  equilibrium  properties  of  gases,  it  is  of  the  utmost  importance  in  the  calcula- 
tion of  inelastic  scattering  cross  sections. 
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The  exp-six  potential  was  therefore  abandoned  and  the  Morse  potential 
(76)  used*  The  dependence  of  the  cross  section  on  the  phase  shift  5  in  (78)  is 
shewn  in  Figures  10  and  11.  Two  representative  values  of  energy  are  shovm,  and 
two  values  of  the  angular  momentum  at  each  energy.  These  curves  make  clear  that 
we  have  to  integrate  the  radial  equation  very  acciirately-  in  order  to  obtain  the 
correct  value  of  the  cross  section.  An  inaccuracy  of  Just  1  in  the  phase  shift 
of  the  final  radial  function  can  cause  an  error  of  as  much  as  lO/o  in  the  calculated 
cross  section. 

Finally,  we  have  taken  two  different  values  of  the  parameters  a  and  d 
in  (76)  and  calculated  the  cross  section  as  a  function  of  energy  for  each  of  these 
values.  These  results  are  shown  in  Figures  12  and  13.  The  cross  section  for  the 
base  potential  (76)  is  shown  for  comparison.  The  variation  in  the  phase  shift 
of  the  radial  functions  with  these  changes  in  the  parameters  is  of  the  order  of 
2-3  radians.  This  accounts  for  the  large  changes  in  the  cross  section  with  re- 
latively small  changes  in  the  parameters. 

It  is  seen  by  comparing  Figures  12  and  13  with  Figure  9  that,  as  one 
would  expect,  the  inelastic  cross  section  is  more  sensitive  to  variations  of 
the  parameters  in  the  perturbation  potential  than  to  the  same  variations  in  the 
unperturbed  potential,  A  change  in  the  steepness  parameter   a  of  about  $/o   in 
the  perturbation  potential  prodnces  the  same  effect  as  a  change  of  about  10 /o  in 
a  in  the  unperturbed  potential. 

3.3  Summary 

We  have  seen  that  the  calculation  of  the  inelastic  scattering  cross  section 
requires  the  solution  of  the  radial  wave  equation  for  very  many  values  of  the  angu- 
lar momentum  Si  .  Tnls   requirement  is  related  to  the  large  wave  number  k  (or  small 
de  Broglie  wavelength)  of  the  system.  In  general,  values  of  angular  momentiim  up  to 
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Inelastic  cross-section  as  a  function   of  phase 
stiift  in  final  radial  function. 
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Inelastic  cross -section  as  a  function  of  phase 
shift  In  final  radial  function. 
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about  kd  where  d  is  the  range  of  the  potential,  must  be  considered.  The  radial 
equation  can  not  generally  be  solved  in  closed  form  for  all  $,    except  for  the 
simplest  and  usually  physically  unrealistic  potentials*  Thus  one  is  forced  to  do 
laborious  computations,  ^"his  presents  no  special  difficulty  in  principle. 

An  analytic  solution  to  the  radial  equation  is  available  for  the  s-wave, 

A  =  0,  if  the  potential  is  not  too  complicated.  This  solution  is  the  one  that 
has  been  \ised  until  recently  for  molecular  scattering  problems.  Zener  and  Devon- 
shire have  given  the  transition  probability  for  the  case  Jj  =  0  with  a  Morse  po- 
tential. ¥e  have  compared  the  present  results  with  those  of  Zener  in  Section  2, 

§  2,3  at  one  very  low  energy  (E.'^  7.10   ev) .   It  woixLd  be  of  interest  to  extend 
that  calculation  and  compare  the  results  over  a  larger  range  of  final  energies. 

We  hiave  shown  that  an  accurate  calculation  of  the  inelastic  scattering 
cross  section  requires  a  mogt  precise  knowledge  of  the  details  of  the  interaction 
of  the  molecular  system.  The  value  of  the  inelastic  cross  section  is  much  more 
sensitive  to  the  details  of  the  perturbation  potential  than  to  the  unperturbed 
potential.  This  sensitivity  of  the  inelastic  cross  section  to  the  parameters  of 
the  potentials  is  due  to  the  small  de  Broglie  wavelength  of  the  system.  Because 
of  the  large  reduced  masses,  the  wavelength  is  very  much  smaller  than  the  range  of 
the  potentials  at  the  energies  of  interest  (0«5  -  2  ev).  Consequently  the  transi- 
tion probability  depends  critically  on  the  potentials. 

The  difficulty  i-iith  the  sensitivity  of  the  inelastic  cross  section  to 
the  precise  shape  of  the  perturbation  potential  is  not  easily  resolved.  The  per- 
turbation potential  is  not  'welL  defined  and  can  be  obtained  from  an  expansion  of  the 
unperturbed  potential  only  if  the  molecules  do  not  approach  too  closely.  Inasmuch 
as  a  pertiirbation  scheme  is  absolutely  necessary  for  this  many  body  problem,  it  is 
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clear  that  nore  work  must  be  done  to  find  the  correct  correlation  of  the  internal 
and  external  motions  of  the  molecules. 

VJe  have  also  seen  in  Section  1  that  the  details  of  the  unpertiirbed  poten- 
tial for  most  molecular  systems  of  interest  are  not  known  with  arycertainty.  Most 
of  the  work  on  the  empirical  fitting  of  a  potential  form  to  experimental  data  has 
concentrated  on  the  region  near  V  =  0.  However^no  work  has  yet   been  done  that 
gives  the  repulsive  wall  of  the  potential  acctirately.  It  is  just  this  part  of 
the  potential  that  is  needed  for  the  calculation  of  an  inelastic  scattering  cross 
section.  Furthermore,  since  the  potential  is  not  known  acc\irately  in  this  region, 
the  present  calculations  can  be  believed  only  to  witliin  an  order  of  magnitude. 
The  relative  variation  of  the  cross  section  with  the  change  in  the  potentials  is 
correct;  however,  the  actual  values  for  the  cross  section  are  not  accurate. 

Here  we  note  that  the  calculation  of  the  inelastic  scattering  cross  sec- 
tion can  itself  be  used  to  discover  the  form  of  the  intermolecular  potential  in 
the  region  from  r  =  0  to  the  zero  of  the  potential.  Since  the  value  of  the  cross 
section  is  acutely  sensitive  to  the  precise  shape  of  the  potential  in  this  region, 
one  can  hope  to  fit  data  on  inelastic  scattering  very  precisely  by  adjusting  the 
parameters  of  the  potential.  Unfortunately,  very  scant  data  are  available  on 
molec\ilar  scattering,  >diether  elastic  or  inelastic.  Such  a  program  must  await  the 
availability  of  experimental  molecular  scattering  data. 
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